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ABSTRACT 

Background: MicroRNAs (miRNAs) are small non-coding RNA molecules that control the post-transcriptional gene 
expression. They play a pivotal role in the regulation of important physiological processes. Variations in miRNA 
genes coding for mature miRNA sequences have been implicated in several diseases, however, its role in type 1 
diabetes mellitus (TIDM) is still underestimated. Objective: To investigate the association of miR-196a2 
(rs11614913) and miR-146a (rs2910164) polymorphisms with T1DM in Egyptian children and adolescents. 

Patients and Methods: A case-control study was done on 100 Egyptian children diagnosed with T1DM and 
additionally, 100 healthy age- and sex-matched subjects non-diabetic (control group) in the Department of Pediatrics 
and Endocrine and Metabolism Unit, Zagazig University Children Hospital from March 2021 to August 2021. 
Random blood sugar (RBS), glycosylated hemoglobin (HbA 1c) and laboratory assessment for lipid profile were done 
for all patients. 

Results: There was no significant difference between CC, CT and TT genotype in age, BMI, sex, onset, and disease 
duration, HbAlc, TC, LDL or HDL. There was a significant increasing in TG level in TT genotype compared to other 
genotypes. There was no significant differences between different genotypes in any of the studied parameters. 
Conclusion: The present study demonstrated that variants rs11614913 T/C and rs2910164 G/C were linked with the 
risk of TIDM. The data suggested that rs11614913 T/C and rs2910164 G/C could be considered as novel risk factors 
in the pathogenesis of T1DM in the Egyptian population. 
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INTRODUCTION upregulated 9-fold in cells transfected with miR- 

Type 1 diabetes mellitus (T1DM) results from 196a2-C allele but increased only by 4.5-fold with 
a cellular mediated autoimmune destruction of miR-196a2-T allele. Other reporter demonstrated that 
pancreatic B-cells, leaving patients insulin-dependent increased expression of miR-196a2 was associated 
for life. The triggering of autoimmunity against B-cells with the C allele ©. 
arises from a multifaceted interaction between multiple Currently, different case control studies have 
genetic and environmental risk factors ®. put the spotlight on miR-196a2 also in various human 

miRNAs are key regulators of gene expression diseases regarding its role in developmental control, 
and, like other genes, their coding sequences are regulation of immune response and regulation of cell 
subject to genetic variation. Polymorphisms in miRNA differentiation and migration. Jameel et al. © study 
genes can have marked effects on miRNA found that miR-196a2 polymorphism was associated 
functionality at all levels, including miRNA with diabetes disease. The study of Assmann et al. ” 
transcription, maturation and target specificity and, indicates that miR-146a rs2910164 polymorphisms are 
consequently, they may contribute to TIDM associated with protection for TIDM in Brazilian 
pathogenesis ®. Recent studies have reported that populations. These associations are biological 
miRNAs are critical posttranscriptional regulators of plausible considering the involvement of miR-146a in 
30-50% of human gene expression and thus they may immunity and inflammation, which are key players in 
contribute in different pathological conditions. The T1DM pathogenesis. 
SNPs in the sequences of miRNAs genes are relatively The aim of this study was to investigate the 
rare and highly conserved, which indicates the association of miR-196a2 (rs11614913) and miR-146a 
functional importance of these SNPs. The miR-196a2 (rs2910164) polymorphisms with TI1DM in Egyptian 
rs11614913 C > T polymorphism is located on children and adolescents. 
chromosome 12 (12q13.13) ®. 

To date, miR-196a gene polymorphism of PATIENTS AND METHODS 
rs11614913 is one of the most investigated SNPs in A case-control study design was conducted on 
case control studies concerning cancer and cardiac 100 Egyptian children diagnosed with TIDM aged 
diseases but no previous study has reported its from 3 to 18 years and 100 healthy age- and sex- 
association with TIDM ™. The association of miR- matched subjects were selected as non-diabetic group 
196a rs11614913 with the previous diseases could be (controls) recruited from the Department of Pediatrics 
explained by findings reported by Hoffman and his and Endocrine and Metabolism Unit, Zagazig 


colleagues © that miR-196a2 expression was 
© O This article is an open access article distributed under the terms and conditions of the Creative 
© Commons Attribution (CC BY-SA) license (http://creativecommons.org/licenses/by/4.0/) 
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University Children Hospital from March 2021 to Aug 
2021. 


Ethical considerations: 

Written informed consent was obtained 
from all children’s’ parents or their relatives and 
the study was approved by the Research Ethical 
Committee of Faculty of Medicine, Zagazig 
University (International review board IRB#:6632- 
3-1-2021). The study was done according to The 
Code of Ethics of the World Medical Association 
(Declaration of Helsinki) for studies involving 
humans. 


Inclusion criteria: Patients diagnosed with T1DM 
according to the criteria of American Diabetes 
Association guideline, and all included cases were 
diagnosed as TIDM using C peptide protein test and 
were under regular treatment with human insulin 
therapy. 

Exclusion criteria: T2DM or specific diabetes. 
Detailed history was taken from included patients 
including the following: 

Full history taking included: Age and sex, detailed 
medical history including; chronological age, duration 
and age at onset of diabetes, type and dose of insulin, 
compliance to diet and insulin therapy, detailed dietetic 
history in addition to family history of diabetes, 
hypertension, coronary heart disease and stroke, and 
history suggestive of diabetic microvascular and 
macrovascular complications. 

Thorough clinical examination including height, 
weight, waist circumference (WC), body mass index 
(BMI=weight in kg/height? in meters). 


Biochemical assessment: Random blood sugar (RBS), 
glycosylated hemoglobin (HbAlIc), Laboratory 
assessment for lipid profile including: serum total 
cholesterol (TC), serum triglycerides (TG), high 
density lipoprotein-cholesterol (HDL-C) and low 
density lipoprotein-cholesterol (LDL-C) = TC-(HDL- 
C)-(TG/5). 


Table (1): Demographic data of the studied groups 
Variable 
Age (years) 


Mean+SD 
Range 


12.55 + 3.56 
3 - 18 


Control 
(n=100) 


13.19 + 2.77 
4-17 


Genotyping and analysis of polymorphisms: 
Expression of miR-146a: 

Plates were positioned in a real-time PCR 
thermal cycler (ViiA7 Real-Time PCR System; 
Thermo Fisher Scientific (USA) and heated for 10 min 
at 95°C, followed by 50 cycles of 95°C for 15 s and 
62°C for 1 min. Only polymorphisms that have minor 
allele frequencies higher than 1% was considered. 


Expression of miR-196a2: 

RNA is isolated using RNeasy Mini Kit 
(Qiagen, Hilden, Germany) according to the 
manufacturer’s instructions, then the extracted miRNA 
was reverse transcribed using a miScript reverse 
transcription (RT) kit (Qiagen, Hilden, Germany). 
Storage of cDNA specimens was done at —80 °C. 
Further qRT-PCR analysis was performed using 
miScript SYBR green PCR kit (Qiagen, Hilden, 
Germany). For the quantification of miR-196a2, we 
used Hs_miR-196a*_1 miScript Primer Assay (Assay 
ID MS00008960, Qiagen, Germany). QPCR was run 
on Applied Technologies’ Stratagene Mx3000P 
(USA), according to the manufacturer’s 
recommendation. 


Statistical Analysis: 

The collected data were entered to and 
analyzed by computer using Statistical Package of the 
Social Sciences version 25 (SPSS). Results were 
presented by tables and graphs. Quantitative data were 
presented as mean and standard deviation and were 
compared by independent t-test or one way ANOVA 
test. Qualitative data were presented as frequencies and 
proportions. Pearson Chi square test (y°) or Fisher’s 
exact were used to analyze qualitative independent 
data. P value of <0.05 was taken as significant. 


RESULTS 

Table 1 shows that there was no statistical 
significance difference between the studied groups as 
regard age, BMI and sex distribution. 


Sig. test P value 





BMI (Kg/m7’) 
Mean+SD 
Range 


17.8 + 4.63 
12.7 — 25.2 


18.3 + 4.90 
10.9 - 25 





Sex (%) 
Male 
Female 





57 (57.0%) 
43 (43.0%) 


SD: Standard deviation t: Independent t test 








48 (48.0%) 
52 (52%) 


x: Chi square test 











Table 2 shows that, the age of onset and duration of TIDM among cases group. 
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Table (2): History of the cases group (n = 100) 


Variable Cases (n = 100) 


Age of onset (years): 


Mean+SD 6.45 + 2.01 


Range 2.5-13 





Disease Duration 
(years): 
Mean+SD 


Range 
SD: Standard deviation 


7.66 + 2.5 
0.5 - 14 








This table shows that, there was a statistical significant 
increasing in HbA1c, cholesterol, TG and LDL among 
cases compared to control group. But there was no 
statistical significant difference between the studied 
groups in HDL level. 





Table (4): Relation between miR-196a2 and miR-146a 
among cases group (n = 100) 


19 
(32.8%) 
16 
(47.1%) 
1 
(12.5%) 

7C:Chi square test **: Highly significant 
Table 5 shows that there was no statistical 
significant difference between CC, CT and TT 


genotype in age, BMI, sex, onset, disease duration, 
HbAIc, TC, LDL or HDL; but there was a statistical 


0.001" 

















significant increasing in TG level in TT genotype 


Table (3): Laboratory findings among the studied compared to other genotypes. 


groups 
Table (5): Relation between miR-196a2 (rs11614913) 


Cases [Control polymorphism and different parameters among cases group 


Variable (n=100) | (n=100) 


P value 
HbAIc Variable 
(%) : <0.001** 
Mean+SD 


Cl 


(mg/dl) 
Mean+SD 


11.964 12.443. |12.84+ 
3.5 18 3.41 








18.05+ | 17.814) 18.34 
5.23 4.9 5.5 








23 19 15 
(40.3%) | (33.3%) | (26.3%) 
17 17 9 
(39.5%) |(39.5%) (20.9%) 








Mean+SD 





5.234 | 5.54 | 6.14 


HDL 1.67 | 181 | 2.02 


(mg/dl) 
Mean+SD 




















6.95+ | 7.314 | 7.84 


2.37 2.25 2.56 





Cl: Cholesterol, TG: Triglyceride, LDL: Low 
density lipoprotein, HDL: High density lipoprotein, 77+ 7.9+ |8.41+ 
SD: Standard deviation, *: Significant, **: Highly Mean+SD | 1.35 1.46 2.8 


significant Cl (mg/dl) | 151.23 |147.23 | 150.24 
Mean+SD | +9.24 | +6.28 | +34.6 
TG (mg/dl) 82.334 |83.83+ |98.16+ 
Mean+SD | 6.97 5.37 30.01 
LDL 
(mg/dl) | 100.34 | 105+ |103.44 
Mean+SD | 3.45 3.19 +3.5 
HDL 
(mg/dl) 51.94 | 49.14 | 50.01 
Mean+SD | 8.16 | +7.96 | +7.94 


Test: F ANOVA test and xê: Chi square test *: Significant. 








Table 4 shows that there was a statistical 
significant increase in frequency of CC genotype of 
miR-146a among cases that had CC genotype in miR- 
196a2. However, there was a statistical significant 
increasing in frequency of GG genotype of miR-146a 
among cases that had TT genotype in miR-196a2. 
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DISCUSSION 

Growing evidence indicates that microRNAs 
(miRNAs) have a key role in processes involved in 
TIDM pathogenesis, including immune system 
functions and B-cell metabolism and death. Although 
dysregulated miRNA profiles have been identified in 
T1DM patients, results are inconclusive; with only few 
miRNAs being consistently dysregulated among 
studies ®. The miR-146a rs2910164 polymorphisms 
were associated with protection for T1DM in Brazilian 
children. These associations are biological plausible 
considering the involvement of miR-146a in immunity 
and inflammation, which are key players in TIDM 
pathogenesis . 

To date, no study has evaluated the association 
between polymorphisms in miR-146a, and miR- 
196a2 genes and T1DM development. So, the current 
study aimed to investigate the association of miR- 
196a2 (rs11614913) and miR-146a (1rs2910164) 
polymorphisms with the susceptibility of T1DM in 
Egyptian children and adolescents at Zagazig 
University Children Hospital. The study population 
consisted of 100 Egyptian children diagnosed with 
T1DM and 100 healthy age- and sex-matched subjects 
were selected as non-diabetic group (controls). 

The current study revealed that, there were no 
statistical significant differences between the studied 
groups in age, BMI or sex distribution. Mean age and 
mean BMI were insignificantly higher in control 
subjects than in T1DM patients. These results were in 
accordance with the study of Assmann et al. . 

However, the age of onset of TIDM among 
cases group ranged from 2.5 to 13 years with mean 
6.45 years while duration of diabetes ranged from 6 
months to 14 years with mean 7.66 years. Similarly, 
Hassan et al. ® reported that the mean age at onset of 
diabetes in the studied group was 8.2+ 2.6 (range 2.5— 
14.5) years, and the mean duration of diabetes was 4.3 
+2.7 (range 1.0—12.0) years. 

Regarding the laboratory findings among the 
studied groups, the current study revealed that, there 
was a Statistical significance increase in HbA\Ic, 
cholesterol, TG and LDL among cases with T1DM 
compared to control group; but there was no statistical 
significance difference in HDL levels between the 
studied groups. These findings were in agreement with 
Hassan et al. ® who reported that, dyslipidemia was 
found in a significantly higher percentage of children 
and adolescents with T1DM (65%) compared to that of 
the non-diabetic control group (28.2%). 

Regarding the relation between miR-196a2 
and miR146a genotype polymorphism, the current 
study revealed that, there was a statistical significant 
increase in frequency of CC genotype of miR-146a 
among cases that had CC genotype in miR-196a2. 
Also there was a statistical significand increase in 
frequency of GG genotype of miR-146a among cases 
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that had TT genotype in miR-196a2. These findings 
were going with the study of Khan et al. %. 

A previous study has shown that miR-196a 
rs11614913 T/C is linked with DM, Downregulation of 
its mature miRNA may play a pivotal role in TI1DM. 
The study of Khan et al. “ found that the genotype 
TT (n = 76/22.4%) of rs11614913 T/C was more 
frequent in cases as compared to controls (n = 
33/13.98%). On the other hand, the genotype CC was 
more common in the control group (n = 130/38.46%) 
than patients (n = 84/35.59%). In patients, 
the heterozygote CT was more common (n = 
178/52.66%) than the control group (73/21.59%). 

Comparison of alleles C and T showed a 
reduced risk for the C allele (OR = 0.4377). Following 
the results of Ibrahim et al. ©, in a study carried out 
in Cairo, Egypt showed that the frequency of TT 
genotype was more common in patients than the 
healthy individuals. In support of this study, this study 
also found a significantly higher frequency of the T 
allele in patients than in controls (41.7% vs. 31.3%), 
while the C allele was more frequent in controls. The 
study stated that the association of the TT genotype 
with TIDM remained significant. The study also 
reported decreased relative expression of miR-196a2 in 
patients compared to controls. 

A previous study suggested that miR-146a 
plays a crucial role in the prognosis of DM by 
contributing to the metabolism, proliferation, and death 
of B-cell. It can be detected in serum, T cells, and a- 
and P-cells of patients with DM, supposing that the 
miR-146a is a therapeutic target and potential 
biomarker “”, A previous study demonstrated that 
rs2910164 G/C leads to an unstable pre-miR-146a 
structure °), The study of Khan et al. was carried 
out in 94 patients and 69 healthy individuals. By 
applying a co-dominant statistical model GG vs GC vs 
CC, a significant difference in genotype frequencies 
was observed (p-value = 0.0001). The genotype GG 
was more frequent in control group (62.3%) as 
compared to diabetes patients (13.8%), whereas the 
genotype GC was more common in patients (52.1%) as 
compared to control (23.1%). The homozygote CC 
was also more common in patients (34.0%) than in 
control (14.5%). The genotype GG reduced the risk of 
diabetes with an odds ratio of 0.2439 (CI = 0.1524— 
0.3902), whereas the recessive model showed an 
increased risk for the genotype CC OR = 2.123 (CI = 
1.118—4.033). 

Khan et al.’s report regarding the 
comparison of the two alleles (G and C) showed a 
reduced risk for G allele 0.2439 (CI = 0.1524—0.3902). 
Thus, it can beconcluded that the C allele of 
rs2910164 G/C in miR-146a increases the risk of DM. 

Regarding the relation between miR196a2 and 
miR146a genotypes and different parameters 
(demographic characteristic data and laboratory 
findings), the current study revealed that, there were no 
statistical significant differences between CC, CT and 
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TT genotype in age, BMI, sex, onset, duration, HbAIc, 
TC, LDL or HDL but there was a statistical significant 
increase in TG level in TT genotype compare to other 
genotypes. 

Age, gender, age at TIDM diagnosis, BMI, 
lipid profile, and HbAIc were similar between subjects 
carrying the minor alleles of the miR- 
146a rs2910164, and miR-196a2 rs11614913 
polymorphisms, analyzed individually, and subjects 
carrying ancestral homozygous genotypes (P > 0.05). 
However, TIDM patients carrying the miR- 
196 rs11614913 C allele (dominant model) had lower 
fasting plasma glucose levels compared to patients 
with the G/G genotype (127.0 + 65.4 vs 
180.4 + 107.6; P=0.001). In non-diabetic subjects, 
age, gender, HbAlc, and BMI, were similar between 
different genotypes of the analyzed polymorphisms 
(dominant model, all P > 0.05) ™. 


CONCLUSION 

The present study demonstrated that variants 
rs11614913 T/C and rs2910164 G/C were linked with 
the risk of TIDM. The data suggested that rs11614913 
T/C and rs2910164 G/C could be considered as novel 
risk factors in the pathogenesis of TIDM in the 
Egyptian population. 


RECOMMENDATIONS 
Further studies performed on larger sample sizes 
and another independent Egyptian sample, are 


recommended to validate and confirm our findings 
among Egyptian population. Moreover, additional 
studies are needed to clarify the functional effect of the 
rs2910164 and rs11614913 polymorphism in immune 
cells as well as pancreatic B-cells. 
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